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Abstract 

We derive several relations between CP violating rate differences A(i? — > 

PV) = T{B PV) - T{B PV) for charmless B ^ PV decays in the 

Standard Model using SU(3) flavor symmetry. It is found that although 

the relations between branching ratios of AS* = and AS* = — 1 processes 

are complicated, there are simple relations independent of hadronic models 

between some of the A5 = and A5 = — 1 rate differences due to the uni- 

tarity property of the Kobayashi-Maskawa matrix, such as A(S ■n'^ p~) = 

-A{B ir+K*-), A{B 7r-p+) = -A{K-p+). SU(3) breaking effects 

are also estimated using factorization approximation. These relations can be 

tested at B factories in the near future. 
PACS numbers: 13.20.He, ll.30.Er, 11.30.Hv 
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I. INTRODUCTION 



Several charmless two body decay modes of B^^d niesons have been observed at CLEO 
|l|J^. These data have provided interesting information about the Standard Model (SM) 
0-0. With increased luminosities for B-factories at CLEO, KEK and SLAC, more useful 
information about charmless Bu ^. decays will be obtained. The SM can then be tested in 
more details. At present only some weak upper limits on the branching ratios have been 
obtained |^ for charmless Bg decays. However, more data on Bg decays will become available 
from hadron colliders, such as CDF, DO, HERAb, BTeV and LHCb in the future. These data 
will help to further test the SM [§. Theoretical predictions, on the other hand, are limited 
by our inability to reliably calculate hadronic matrix elements related to B decays although 



progresses have been made in recent years [|T0|- In the lack of reliable calculations, attempts 
have been made to extract useful information from symmetry considerations. SU (3) flavor 



symmetry |pT|-[r^ is one of the symmetries which has attracted a lot of attentions recently. 
It has been shown that using SU(3) symmetry it is possible to constrain and to 

determine 0] one of the fundamental parameters 7 in the SM for CP violation by measuring 
several charmless hadronic B meson decay modes. SU(3) flavor symmetry also predicts many 
interesting relations between CP violating observables in the SM, such as rate differences 
between AS = and = — 1 B decays. These relations will provide important test for 
the SM. 

CP violating rate differences A{B 0102) = r(i? 0102) — r(i? 0102) for B decays 
have been studied before |T5Hl7| . Here Oj can be one of the octet pseudo-scalars P or one of 
the octet vector mesons V, respectively. A complete study for the charmless B PP case 
was carried out in Ref. [T^. For B —>■ VV, as long as rate differences are concerned, the 
situation is similar to B ^ PP. For B PV, the situation is more complicated because 
there are more independent SU(3) invariant decay amplitudes and more relations exist for 
these decays. In this paper we will concentrate on CP violating rate difference relations for 
B PV decays in the SM using SU(3) flavor symmetry. B PV decays using SU(3) 
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flavor symmetry have also been studied in the literature recently in Ref. p8|JT9|] with different 
emphases. Keeping the leading order contributions, Ref. |]T8[ studied available data from 
CLEO and obtained information about the phase angle 7 and proposed new tests for the 
SM. Keeping all contributions Ref. studied electroweak penguin effects and information 
about 7. In our analysis we also keep all contributions, but with the emphasis on possible 
relations between rate differences which have not been studied before in the literature for 
B to PV decays. At present CLEO collaboration has measured several B PV modes. In 
the near future, B factories will measure many more of these decays with more precision and 
rate differences may be measured. CP violation in the SM will be tested using B — > PV 
decays. 

The paper is arranged as following. In section II we present the SU(3) decay amplitudes 
for charmless B — > PV decays. In section III we study relations between rate differences for 
various B to PV decays. In section IV, we estimate SU(3) breaking effects on the relations 
obtained. And in section V we draw our conclusions. 

II. 5[7(3) DECAY AMPLITUDES FOR CHARMLESS B ^ PV 

The quark level effective Hamiltonian up to one loop level in electroweak interaction for 
charmless hadronic B decays, including the QCD corrections to the matrix elements, can be 
written as 

Hi^f = ^[k,k;(ciOi + - + Vt,v,icrm (i) 

i=3 

The coefficients Ci 2 and cl^ = cl — c^, with j indicates the internal quark, are the Wilson 



Coefficients (WC). These WC's have been evaluated by several groups |2^, with |ci,2| >> 
|c^|. In the above the factor VcbV*^ has been eliminated using the unitarity property of the 



KM matrix. The operators Oi are defined as [20 
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Oi = {qiUj)v~A{uihj)v-A , O2 = {qu)v~A{ub)v-A , 

03,5 = {qb)v^A Eg'{q'q')vTA , 04,6 = {qibj)v-A T.q'{q'jq'i)vTA , 

^7,9 = l{qh)v-AY.q' (iq'{q'q')v±A , Os.io = |(gi&i)y-AEg'eg'(g^g-)y±A , 

On = + 75)& , Ou = ^qa,,F'^''{l + ^,)h. 

where {ah)v-A = a7^j(l — 75)^, G^'^ and F^'^ are the field strengths of the gluon and photon, 
respectively. 

At the hadronic level, the decay amplitude for B PV can be generically written as 

AiB PV) =< P V\HyB >= V^,V:^T{q) + Vt,V,\P{q) , (3) 

where T(g) contains contributions from the tree operators Oi,2 as well as penguin operators 
03_i2 due to charm and up quark loop corrections to the matrix elements, while P{q) 
contains contributions purely from penguin due to top and charm quarks in loops. We 
would like to clarify the notation used here. The amplitude T in eq. ^ is usually called the 
"tree" amplitude which will also be referred to later on in the paper. One should, however, 
keep in mind that it contains the usual tree current-current contributions proportional to Ci,2 
and also the u and c penguin contributions proportional to c" — with i = 3 — 12 which is 
small due to cancellation between contributions from up and charm quarks in loops. Also, in 
general, it contains long distance contributions corresponding to internal u and c generated 
intermediate hadron states |2T|. In our later analysis, we do not distinguish between the 



tree and the penguin contributions in the amplitude T unless specifically indicated. 

The relative strength of the amplitudes T and P is predominantly determined by their 
corresponding WC's in the effective Hamiltonian. For AS = charmless decays, the dom- 
inant contributions are due to the tree operators Oi^2 and the penguin operators are sup- 
pressed by smaller WC's. Whereas for AS* = — 1 decays, because the penguin contributions 
are enhanced by a factor of VtbVf*/VubV*g ~ 55 |^ compared with the tree contributions, 
penguin effects dominate the decay amplitudes. In this case the electroweak penguins can 



also play a very important role p2|. In all AS* = and AS* = — 1 processes, both the tree 
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and penguin contributions have to be present in order to have rate differences. One should 
carefully keep track of all different contributions. 

The operators Oi,2, Os-e.n, and 07-10,12 transform under SU(3) symmetry as 3a + 3;, + 
6 + 15, 3, and 3a + 3^ + 6 + 15, respectively. These properties enable us to write the decay 
amplitudes for B PV in only a few SU(3) invariant amplitudes. 

For the T{q) amplitude, for example, we have [|l^ 



T{q) = A-,{T)B,H{l)\VtM[) 

+ Cy{T)BiVlM^H{2,y + C^{T)BiMiV^H{?,y 

+ Al{T)B,H{QriV^^Mf + A^'{T)B,H{QyiM^^/t 

+ C^iT)B,V;Hi6)i'Mi + C','{T)BM]H{QtVl 

+ AUT)B,H{T^tV^^Mf + A^^{T)B,H{T^riM^^/t 

+ CUT)B.V;Hm'Mi + C§iT)B,MiHm'vl , (4) 

where Bi = Ba, B^) = {B-,B^, 5°) is a SU(3) triplet, M/ and are the SU(3) pseudo- 
scalar and vector meson octets, respectively. In our analysis, we will also include the SU(3) 
singlets for both the psuedo-scalar and vector mesons, such that M and V become the U(3) 
nonet, to have some idea about the decay amplitudes for B decays involving these particles. 
The SU(3) invariant amplitudes Ai{T) and Ci{T) are in general complex due to final state 
interactions. The matrices H{i) contain information about the transformation properties of 
the operators Oi_i2. 

For q = d, the non-zero entries of the matrices H[i) are given by [jl2 



H(15)l' = H(15)l' = 3 , HiT5)f = -2 , H(T5)f = H(15)f = -1 . (5) 
And for g = s, the non-zero entries are |15] 



H(15)l' = H(15)l' = 3 , Hil5)f = -2 , H(15)f = H{\5)f = -1 . (6) 
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There are similar amplitudes for the penguins. We will indicate these SU(3) invariant 
amplitudes by Ai{P) and Ci{P). 

The decay amplitudes can be written in terms of the SU(3) invariant amplitudes as 

AiB -> PV) = Vu,V:XAnT) + a^AfiT) + h't Cj (T) + 6f (T)) 

+ V,,V:^{aUUP) + afAf{P) + h'^ Cr {P) + &f C^f (P))- (7) 

The index i is summed over i = 3, 6, 15. For i = 3, there is only one amplitude. We use 
the convention: A^ = A^, = a^, and a^^A^^ = 0. 

One can easily obtain the decay amplitudes for B PP and B —>■ VV from the above 
by replacing VJ by Mj, and Mj by VJ in the above expression, respectively. The SU(3) 
invariant amplitudes are then replaced by Ai = AY + Af and Ci = CY + C*^. In these 
two cases, due to the anti-symmetric nature in exchanging the upper two indices of H]^{Q) 
and the symmetric structure of the two mesons in the final states, Cq — A^ always appear 



together ||T2[ if the singlets in the final states are removed. Therefore in these cases there 
are in total 5 independent SU(3) invariant amplitudes for each case. However if singlets 
are included, decay modes involving singlets in the final states can be used to separate Aq 
and Cq. In this case there are in total 6 independent amplitudes. For B PV, even the 
symmetric structure of the two mesons in the final states is lost, there are relations between 



y4g and Cg if the singlets r/i and the combination of {J2/3)uj + (p/V^ are absent, there are 



total 10 independent invariant amplitudes which agree with the analysis in Ref. [|T^. When 
the singlets are introduced, A^ and Ci are all independent. There are total 11 of them. 
The analysis in this case is more complicated compared with the cases for B PP and 
B VV. Expanding eq. ^, we obtain the coefficients for each individual decays. In tables 
1 - 3 we list the coefficients and bi for AS* = B ^ PV decays, and in Tables 4-6 for 
AS = -1 B ^PV decays. 

Two remarks are in order: a) The amplitudes Ai correspond to annihilation contributions, 
as can be seen from eq. ^ where Bi is contracted with one of the index in H{j), are small 
compared with the amplitudes Ci from model calculations. The smallness of these annihila- 



tion amplitudes can be tested using Bd K~K*~^, K^K*~ and Bg vr+p", it~ p'^, n^p^ , ti^uj 
because these decays have only annihilation contributions, b) Many analyses have been car- 
ried out using SU(3) classification of quark level diagrams ||13[ in the literature. In most 
cases such analyses give the same results as the use of SU(3) invariant amplitudes discussed 
here. However, in some cases the classification according to quark level diagrams without 
care could be misleading. One should be careful to include all possible contributions to 
have a complete study For example, neglecting annihilation contributions, the op- 

erators Oi,2 would have vanishing contributions to Bu — ^ K~ K*^ , K*^ , K*^ , p~ , 
Bd K^K*^,K^K*^, and 5, K^K*^,K^K*^. However, from the Tables, we find that 
these decays are proportional to — C^'^^ — Cjg-*^- These combinations are not neces- 

sarily zero without additional assumptions. Whether naive diagram analysis is valid has to 
be tested experimentally. This can be achieved by measuring the decay modes mentioned 
in the remarks a) and b). Without evidence from experimental data, we have to keep the 
most general form and carry out analyses accordingly. 

III. CP VIOLATING RATE DIFFERENCES FOR CHARMLESS B ^ PV 

One can find many relations between different B decays using SU(3) flavor symmetry to 
test the SM in different ways. A particularly interesting class of relations for B PV is 
the CP violating rate difference A{B PV) = T{B ^ PV) - T{B PV). Theoretical 
calculations of rate differences are difficult because one not only needs to calculate the short 
distance decay amplitudes but also the long distance contributions, especially the final state 
interaction phases. This is a difficult task at present. Due to this reason, even though 
CP violation is measured for certain B decay modes and compared with hadronic model 
calculations, one is not sure if the results agree with the SM predictions. In this section we 
will derive several relations which do not depend on the dynamics of hadronization processes 
using the SU(3) decay amplitudes obtained in the previous section. These relations can 
provide hadronization model independent tests for CP violation in the SM. 
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Using the SU(3) decay amplitudes obtained in the previous section, one can find that 
some decay amphtudes for AS' = and AS = —1, B ^ PV decays have the following 
peculiar form [ll-^ig , 

Aid) = v^,v:,T + Vt,v:,p, 

Ms) = Vu,V:,T + Vt,VZP. (8) 

Due to different KM matrix elements involved in A{d) and A{s), although the amplitudes 
have some similarities, the branching ratios are not simply related. However, when consid- 
ering rate difference, the situation is dramatically different. Because a simple property of 
the KM matrix element 0, /m(KfeKd^i6^td) = in the SU(3) hmit we 

have 

A(d) = -A(s), (9) 
where A(i) = — \A{i)\'^)\ab I {^T^rnB) is the CP violating rate difference defined earlier 



and \ab = Y 1 ~ 2(m^ + mf)/m\ + (m^ — m^)2/m^ with niafi being the masses of the two 
particles in the final state. 

The above non-trivial equality does not depend on the numerical values of the final state 
re-scattering phases. Of course these relations are true only for models with three gener- 
ations. Therefore they provide dynamic model independent tests for the three generation 
Standard Model. 

We find the following equahties: 

(1) A(5„ K'K*"^) = -A(5„ ^ 7r-i^*°) , 

(2) A{Bd = -A{Bs K^K*^) , 

(3) A(5„ K'^K*-) = -A(5„ K'^p-) , 

(4) A{Bd K^K*^) = -A{Bs K^K*°) , 

(5) A(5, ^ 7r^p+) = -A(5, ^ K-K*+) , 

(6) A{Bs ^ 7T-K*-') = -A{Ba ^ p+) , 
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(7) A{Bd- 


TT+p ) = 


A / 7~» 

-A{B, 


K^K* ) , 


(8) A{Bs - 


K+p ) = 


-A{Bd - 




(9) A(S„- 


^ViP ) = - 


A / 7~» 

-A{Bu 


ViK* ) , 


(10) A{Bs - 


- ViK*') = 


-A{Ba - 


^ r/ii^*°) , 


(11) A(S,- 


K-K*+) 


= -HBs 


7r-p+) , 


(12) A{Ba- 


K+K*-) 


= -A{Bs 


-ttV)- 



(10) 



If it turns out that the annihilation contributions arc all small, as can be tested in 
Ba K~K*^, K^K*~ and Bg p~ 7r~p"'", 7r°p°, 7r°a;, there are additional relations for 
rate differences. We find the following equalities, 

(1) = (2) , (3) = (4) , ^^^^ 

(5) = (6) , (7) = (8). 

IV. SU(3) BREAKING EFFECTS 

The relations obtained in the previous section hold in the SU(3) hmit. The SU(3) 
symmetry need to be tested experimentally. To this end we comment that SU(3) pre- 
dictions for B PV with charmed vector meson V can be independently tested by 
using B — > D*T^ and B — > D*K. These decays receive contributions from the tree 
operators Oi^2 only. This provides a clear test for SU(3) symmetry free from penguin 
contaminations. In the SU(3) symmetry, the ratio of these branching ratios is equal to 
r = Br{B D*n)/Br{B D*K) = |Kd/KsP- Factorization calculation gives a SU(3) 
breaking factor and r = (/^//|)|Kjd/KtsP- This can be tested experimentally. Results ob- 
tained from this will give us a good guidance on how well CP violation can be tested in 
charmless B — > PV decays. 

If SU(3) is broken these relations need to be modified. We now study how these relations 
are modified when SU(3) breaking effects are included. Since no reliable calculational tool 
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exists, in the following we will use factorization approximation neglecting the annihilation 
contributions to estimate the SU(3) breaking effects. We find 0,^] 

(2) A(5, K'K*') = -( itZ^I'f J ^A(i?. - K'K*') , 

(3) A(5„ K'Kn = -(S^^!/^)^A(i?„ K'p-) , 

(4) A(5, ^ i^-Oir*") = -{^Z^f^iB, - ^°^*°) ' 

(5) A(5, ^ vr-p-) = -(7^SSt)^A(^^ - ^-^*") > 

JkAq yrriK) 



(6) A(5, ^ 7r-/r+) = -( . \-,.'\ y^iB, ^ ir-p 



(7) A(5, ^ vr+p-) = -( /^^^^^;/; )^A(i?. ^ ir+ir*- 



(8) A(i?. ^ irV) = -( /'^;il/7\ )'A(i?, ^ TT+ir*-). (12) 



In the above we have not listed the corrections for the equalities (9), (10), (11) and (12). 
The corrections to the decay modes in (9) and (10) involving rji are complicated because 
there are two corrections to the amplitudes and also because mixings between rji and r^g. It is 
difficult to use these decay modes to have clear tests. The decay modes in (11) and (12) are 
all pure annihilation type of decays which are zero in the naive factorization approximation 
and likely to be small. We will not discuss these decay modes further here. 

The decay modes in (1), (2), (3) and (4) do not have contributions from O12 in the 
factorization approximation. Therefore their rate differences are expected to be very small. 
However one should be careful in drawing this conclusion, as already mentioned earlier, 
due to the fact that the vanishing contributions from Oi^2 to these modes is not a SU(3) 
prediction even when annihilation contributions are neglected. It may happen that the rate 
differences are larger than expected. We have to wait experimental data to tell us more. 

At e'^e^ B factories the best chances to test the above listed relations are from (5), (6), 
(7) and (8). To a good approximation, we have 
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(6) A{B, vrV) ^ -( ;'Jb^.)J ^ n^K*-) . (13) 



The ratios A^^PimD/AQ^^imj^) and Ff-"^(m2)/Ff ^'^(mf^,) are all close to one[3]. The 
corrections are dominated by the decay constants /j. The relation (a) has, then, a larger 
correction than (b) because fp/ fx* ~ 1 and fn/ fx ~ 0.82. 

The decay branching ratios for Bd 7i^K*~ and the combined Bd — > vr^p^, vr^p^ 
have been measured at CLEO with Br{Bd n+K*~) = (22^^^^) x 10-^ and Br{Bd 
TT+p^ +7r^p+) = (35l|o±5) X 10~^, respectively. There is only an upper bound for Br{Bd — 
K^p^) < 25 X 10~^ at the 90% c.l.. If we normalize the rate differences by a common 
branching ratio, for example, Br{Bd tt^K*~), the normalized rate differences for Bd 
7r^p+, K"p^ are of order a few percent while for Bd 7r+p^, n^K*^ are of order 20% 
P). These can be measured at B factories with good precision. The relations predicted can 
be tested. Even earlier a test can be performed with the combined data 

A{Bd - 7r"p+) + A{Bd - TT+p-) ^ -{^fA{Bd - K~ p+) - (-^)'A(5, ^ n+K*-) . 

JK JK* 

When Bs decay modes are measured, such as at CDF, DO, BTeV, HERAb, and LHCb 
all modes in (5), (6), (7) and (8) can provide good tests for the Standard Model. 



V. CONCLUSIONS AND DISCUSSIONS 

We would like to point out that the relations obtained for the CP violating rate differences 
also hold for any three generation models in which flavor changing and CP violating source 
is solely due to KM matrix, such as multi-Higgs doublet models with neutral flavor current 
conservation by exchanging of Higgs bosons at tree level and model with anomalous three 
gauge boson couplings as new physics source. 

To conclude we find that although the relations between branching ratios of AS* = 
and AS" = — 1 processes are complicated, there are simple relations between some of the 
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AS — and AS — —1 rate differences due to tfie unitarity property of the KM matrix, such 
as A{B TT+p-) = -A{B tt+K*'), A{B 7r-p+) = -A(if-p+). We also estimated 
SU(3) breaking effects using factorization approximation. The relations involving Bd and 
Bu decays can soon be tested by CLEO, Babar and Belle collaborations. While the relations 
involving Bg decays can be carried out at hadron coUiders, such as CDF, DO, BTeV, HERAb 
and LHCb. CP violation can be tested in a model independent way for the Standard Model 
using the relations derived in this paper soon. We emphasis that the importance of these 
tests is their independence from hadronic models. We urge our experimental colleagues to 
carry out such analyses. 
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Table 1. SU(3) decay amplitudes for AS = 0, S„ ^ PV decays. 



Decay Mode 


n- hM. 


K-K*^ 
K^K*- 

Bu mp~ 

Bu ^ ttV" 
Bu 7r~0 

Bu TI""*^ 


(01 3 00 10 -1 -1 ) 

(00 1 31 -1 -1 0) 
^(01 1 3 31 1 -3 133) 
^(0113 3110 103) 
^ ( 1-1 3-3-1 1-1 1 5 3 ) 

(00 0000 10 -1 ) 
^ ( -1 1 -3 3 1 -1 1 -1 3 5 ) 
^(0113 3111-13 1) 
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Table 2. SU(3) decay amplitudes for AS = 0, ^ PV decays. 



Decay Mode 



n- hM. 



5, 

Bd 
Bd 

Bd 
Bd 
Bd 
Bd 
Bd 
Bd 
Bd 
Bd 
Bd 
Bd 
Bd 
B,i 



K-K*+ 
K+K*- 

r]i(f) 



7:\ 

^0 



1 

VE 
1 

2a/3 
1 

2x/3 

1 

1 
V6 

1 

1 
2 

1 
V2 

1 

2 



1 1 0-1-2 1 0-1 0-1 



-2-10 10-1 0-1 
3 -1 000000 
3-2101030 
-2 3010103 
-1 3000000 
2 200010 -1 
1111-3-131 
5 5-1-13-1-3 5 
1 1 1 -1 -1 1 -1 
2-1-111110-101 
0-1-1 5 5-1-10-1 5 
2-1-1 1 11 11 1-5-5 
00000000 -1 01 
-1 -1 5 5 -1 -1 -1 1 5 -1 



1 1 

1 -1 1 

1 -1 

1 -1 

1 1 

-2 -2 

2 -1 
-1 



-1 

-2 
-1 
-1 
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Table 3. SU(3) decay amplitudes for AS = 0, ^ PV decays. 



Decay Mode 



n- hK 



Bs 
Bs 
Bs 

Bs 
Bs 
Bs 
Bs 
Bs 



K+p- 
7r-K*+ 



1 

V2 

1 

%/2 



1 

%/2 

1 

V3 

1 
V6 



0-1 0-11 0-1 1-1-1 

01 1 00-10-105 
0-10-10010-101 
0-10-10010103 
00 -1 -1 101030 
00 1 1 -1 -1 050 
0-1-1-1-11 10-1 -1 
2 -1 2 -1 1 -2 -3 2 3 2 
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Table 4. SU(3) decay amplitudes for AS = -1, PV decays. 



Decay Mode 



n- hM. 



By 
Bu 

By 
By 
By 

By 
By 



K-uj 



1 

V6 

1 

%/3 



1 
V2 

1 

%/2 



1 
V2 



01 3 0010 -1 -1 
1-23-6-210163 

113 3 110 10 3 

01 03 001000 -2 
00 1 3101 -2 34 
00103101032 

00 1 310 -1 -1 

01 03 001 -2 143 



Table 5. SU(3) decay amplitudes for AS — —1, B^ — > PV decays. 



Decay Mode 



Bd 
Bd 
Bd 
Bd 
Bd 
Bd 
Bd 
B, 



n+K*- 
K-p+ 



„_ „V „M „y „M I.V uM uV uM i.V iM 

03 Og Og a-^ O3 O3 Og Og 0-^ Oj^ 



15 15 



1 

1 
%/2 



1 

V2 
1 

1 

V6 



10 1-10 1-214 

0-1 0-11 0-1 0-1 2 

-1 -1 1 -2 
0-10-10010103 

00 -1 -1 101030 

010 100-1-2 141 
0-1-1-1-1110-10 -1 
0-1 2-1 2-2 10-1 6-1 
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Table 6. SU(3) decay amplitudes for AS — —1, Bg — > PV decays. 



Decay Mode 



n- 

^3 C^e ^^6 ^15 



M ly uM hV uM lyiM 

"T5 "3 ^^3 "6 



Bs 
Bs 

Bs 
Bs 
Bs 
Bs 
Bs 
Bs 
Bs 
Bs 
Bs 
Bs 
Bs 
Bs 
5, 



K-K*+ 
K+K*- 



7rV° 



7t\ 

^0 



2 
3 

1 
1 

%/3 
1 

V3 

2^ 

3 



V2 






1 


-2 


-1 





1 





-1 





-1 


1 





-1 


-2 


1 





-1 





-1 





1 





3 


-2 


1 





1 





3 





1 


1 


3 


-1 




















1 


-1 


-1 


3 























-1 


-2 


3 





1 





1 





3 








2 


2 


-1 


-1 








3 


2 


1 


-1 


2 


2 











2 





-4 








1 


1 

















-2 








-2 


-2 


1 


1 











-2 








1 


1 

















1 


1 


-1 


2 


2 











-1 





2 








1 


1 






































-1 





2 





1 


-1 


2 


2 
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